function[beta_k,Xi,V_4,REX]=E_V_03(G_4,ND,NO,Hz,i_7);

[Eig_V,Eig_D]=eig(G_4);
Eig_G=diag(Eig_D);
alph=Hz*log(abs(Eig_G));
beta=Hz*atan(imag(Eig_G)./real(Eig_G));
beta_k=sqrt(alph.^2+beta.^2);
Xi=-alph./beta_k;
beta_k=beta_k/(2*pi);

Ind=find((beta_k<10 & beta_k>0.1 & Xi<1.20 & Xi>0.00 & imag(Eig_G)~=0)); 

NN=length(Ind);
Ind=Ind(1:2:NN);
NN=length(Ind);
[B,ind_1]=sort(beta_k(Ind));
V_1=Eig_V(((NO-1)*ND+1):(NO-0)*ND,Ind(ind_1));

V_2=abs(V_1);
V_3=angle(V_1);
for i_0=1:NN
    V_4(:,i_0)=V_2(:,i_0)/max(V_2(:,i_0));
    V_22=V_2(:,i_0);
    k_0=find(V_22==max(V_2(:,i_0)));
    V_5=abs(angle(exp(i*(V_3(:,i_0)-V_3(k_0,i_0)))));

    for i_1=1:ND

        if V_5(i_1)>pi/2
            V_4(i_1,i_0)=V_4(i_1,i_0)*(-1);
        end
    end
end

if NN>0
    REX=[ones(NN,1)*NN beta_k(Ind(ind_1)) Xi(Ind(ind_1)) V_4']';
else
    REX=[NN; zeros(2+ND,1)];
    V_4=0;
end
